Self-trapped bidirectional waveguides in a saturable photorefractive medium 
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We introduce a time-dependent model for the generation of joint solitary waveguides by counter- 
propagating light beams in a photorefractive crystal. Depending on initial conditions, beams form 
stable steady-state structures or display periodic and irregular temporal dynamics. The steady-state 
solutions are non-uniform in the direction of propagation and represent a general class of self-trapped 
waveguides, including counterpropagating spatial vector solitons as a particular case. 
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PACS numbers: 42.65.Tg, 42.65.Jx, 42.65.Sf 

During the last decade spatial screening solitons p] 
have been considered almost exclusively in copropagation 
geometry. Recent progress in generating optical solitons 
consisting of counterpropagating fields by Cohen et al. 
Q has renewed interest in counterpropagating four-wave 
mixing, extensively studied in the past U'4, 5]. However, 
such geometries inphotorefractive (PR) media are prone 
to instabilities @, 0, Q and are often employed for trans- 
verse optical pattern formation Q. In particular, tem- 
poral instabilities were shown to result in self-oscillation, 
chaos, and bistability UB- It is therefore of importance 
to investigate the temporal behavior of counterpropagat- 
ing self-trapped beams in PR crystals with finite response 
time. Furthermore, one may easily envision interest in a 
stable self-adjustable bidirectional connection of two ar- 
rays of beams across a PR crystal. 

In this communication we derive equations for the 
propagation of beams, similar to the bimodal counter- 
propagating solitons in Kerr media and collisions of 
screening PR solitons propagating in opposite directions 
0, together with a time-relaxation procedure for the 
space charge field determining the refractive index modu- 
lation in PR crystals. Dynamical effects are found impor- 
tant for understanding the behavior of counterpropagat- 
ing beams. We display numerically the temporal forma- 
tion of bright spatial screening vector solitons formed by 
counterpropagating beams, and discuss their interactions 
in (1+1) spatial dimensions. Beyond soliton solutions, we 
introduce a more general class of steady-state induced 
waveguides. Additionally, a situation where the interact- 
ing beams do not converge to a stationary structure, but 
rather alternate between different states, is reported. 

We consider two counterpropagating light beams in a 
PR crystal, in the paraxial approximation, under condi- 
tions suitable for the formation of screening solitons. The 
optical field is given as the sum of the counterpropagat- 
ing waves F exp(ikz + iut) + B exp(— ikz + iuot), k being 
the wave vector in the medium, F and B are slowly vary- 
ing envelopes of the beams. The total light intensity I is 
measured in units of the background light intensity, also 
necessary for the generation of solitons. After averaging 
in time on the scale of response time To of the PR crystal, 
the total intensity is given by: 

1 + 1 = (1 + J ) [1 + e {m exp(2ikz) + c.c.} /2] , (1) 



where I Q = \F\ 2 + \B\ 



2FB*/(1+I ) is the modula- 



tion depth, and c.c. stands for complex conjugation. Here 
the parameter e measures the degree of temporal coher- 
ence of the beams related to the crystal relaxation time: 
for e = 0, i.e. when the relative phase of the beams varies 
much faster than Tn, the beams are effectively incoherent 
(see discussion in J2()- In the opposite case e = 1, the in- 
tensity distribution contains an interference term which 
is periodically modulated in the direction of propagation 
z, chosen to be perpendicular to the c-axis of the crystal, 
which is also the x-axis of the coordinate system. Beams 
are polarized in the x direction, and the external elec- 
tric field E e , necessary for the formation of self-trapped 
beams, also points in the x direction. The electric field 
in the crystal couples to the electro-optic tensor, giving 
rise to a change in the index of refraction of the form 
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where n$ is the unperturbed index, r e ff is the effective 
component of the electro-optic tensor (in this case 7*33), 
and E to t is the x-component of the total electric field. It 
consists of the external field and the space-charge field 
E sc generated in the crystal, E to t = E e + E sc . 

The intensity modulates the space charge field, which 
we represent in the normalized form 



E sc / E e — Eq 



— [Ex exp(2ikz) 



(3) 



where Eq is the homogeneous part of the x-component 
of the space charge field, and E 1 (x,z) is the additional 
slowly varying part of the space-charge field proportional 
to e, |<9 z i?i| <C 2fc|i?i|. It is Eq that screens the external 
field, and E\ is the result of the interference pattern along 
the z-direction. It vanishes together with the intensity 
modulation for incoherent beams, i.e. in the limit e = 0. 

In a simplified approach, one assumes a local, isotropic 
approximation to the space charge field, and looks for a 
solution with saturable nonlinearity E sc — E e /(1 + /). 
Substituting Eqs. and © in this expression, neglect- 
ing higher harmonics and terms quadratic in to, we ob- 
tain as a steady-state solution 
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FIG. 1: Counterpropagating dipole-mode vector soliton (a), 
made out of a fundamental beam propagating to the right 
(b), and a coherent dipole beam propagating to the left (c). 
Coupling strength F = 10/3. 



FIG. 2: Bidirectional waveguide, (a) Total intensity distribu- 
tion; (b) right- and (c) left-propagating beams. Layout as in 
Fig- Q parameters e — 1, and T = 5. Initial peak intensities 
I F = I B = 1. 



Temporal evolution of the space charge field is intro- 
duced by assuming relaxation-type dynamics 

Td t Eo + E = -j^j- , (5a) 
J- + -*o 

FTT1 

TdtE\ + Ei = — - — ■ — — , (5b) 
l + Jo 

where the relaxation time of the crystal r is inversely 
proportional to the total intensity r = tq/(1 + I), i.e. 
illuminated regions in the crystal react faster. The as- 
sumed dynamics is that the space-charge field builds up 
towards the steady-state, which depends on the light dis- 
tribution, which in turn is slaved to the slow change of 
the space-charge field. As it will be seen later, this does 
not preclude a more complicated dynamical behavior. 

Selecting synchronous terms in the nonlinear paraxial 
wave equation we obtain the propagation equations 

id z F + d 2 x F = T [E F + Ej.B/2] , (6a) 
-id z B + d%B = T[E B + E$F/2], (6b) 

where the parameter T — (knoXo) r e ffE e , and we 
use the rescaling x — * x/xo, z —> z/Ljj, (F,B) — > 
(F, B) exp(— iTz). Here xq is the typical beam waist 



and Ld = 2hx\ is the diffraction length ^j- Prop- 
agation equations are solved numerically, concurrently 
with the temporal equations. The numerical procedure 
consists of solving Eqs. JSJ for the components of the 
space-charge field with the light fields obtained at every 
step as a guided modes of common induced waveguide. 
This is achieved by an internal relaxation loop, i.e. nested 
within the temporal loop, based on a beam-propagation 
method for the right- and left-propagating components. 
Both loops are iterated until convergence, which however 
is not necessarily reached in the temporal loop. In that 
case a dynamical state is obtained. 

Head-on collision of the beams with initial soliton 
profiles, after temporal relaxation to a steady-state, re- 
sults in the formation of a counterpropagating soliton 
(not shown), similar to the one found in Ref. |2j. One 
can easily generalize this approach, introducing higher- 
order counterpropagating solitons, similar to the multi- 
hump vector solitons in co-propagating geometry (see, 
e.g. Ref. 0). In Fig. ^ we present a particular case of 
a dipole-mode counterpropagating soliton. Dipole beam 
is launched from the right, and a power-matched single 
beam from the left. Such a bimodal counterpropagating 
soliton has been studied in Ref. H and has been found 



FIG. 3: Incoherent (a)-(c) and coherent (d)-(f) interaction 
of two pairs of counterpropagating beams. The initial offset 
is 4ko for in-phase beams propagating to the right in (b) and 
(e), and 2xo for the out-of-phase beams propagating to the 
left in (c) and (f). Parameters and layout as in Fig. 



to be stable in co-propagating geometry The size of 
data windows in all figures is 10 beam diameters trans- 
versely by 2 diffraction lengths longitudinally. 

Shooting initial beams with arbitrary parameters gen- 
erally leads to z-dependent or non-stationary character 
of the beam propagation. In some domain of the initial 
parameters, for example with the relative angle of beam 
scattering 8 close to it and small initial transverse offset, 
our time-relaxation procedure converges to the station- 
ary in time structures, which we denote as steady-state 
self-trapped waveguides. The formation of a single bidi- 
rectional waveguide is shown in Fig. [21 Two coherent 
Gaussian beams are launched at different lateral posi- 
tions perpendicular to the crystal edges, 9 — it. Both 
beams diffract initially, until the space-charge field is de- 
veloped in time to form the waveguide induced by the 
total light intensity, Fig. |21 (a) , and this induced waveg- 
uide traps both beams, Fig.[5](b) and (c). When the ini- 
tial separation is 4 or more beam diameters, the beams 
hardly feel the presence of each other, and focus into indi- 
vidual solitons. For the separation of 2 beam diameters, 
the interaction is strong enough for the beams to form a 
joint waveguiding structure, as it shown in Fig. [21 

We would like to note here, that for a propagation 



FIG. 4: Unstable self-organized beam structure after 116to 
temporal steps, at the moment of symmetry-breaking. Three 
in-phase beams propagate to the right (b), two out-of-phase 
beams to the left (c). The components have equal powers and 
the coupling F = 10. Other parameters as in Fig. Q 

distances exceeding some threshold value, i.e. for larger 
crystal lengths, we observe a longitudinal modulational 
instability developing in time even for the initial beams 
corresponding to the exact steady-state solitons. Mod- 
ulational instability is a topic of ongoing research and 
beyond the frame of present paper. 

Having found different steady-state self-trapped struc- 
tures, we examine the difference between the coherent 
and incoherent interaction of beams. Two steady-state 
solutions with the same boundary conditions but for dif- 
ferent degrees of mutual coherence e are shown in Fig. 01 
Counterpropagating beam components made of two pairs 
of beams are launched with a lateral offset. The beams 
to the right are in-phase, and aim at the center of the 
opposite crystal face. The beams to the left are out-of- 
phase, and launched in parallel. Figs.|31(a)-(c) depict the 
incoherent interaction, e = 0. The beams attract, focus 
and overlap tightly, but the ones to the right (b) are still 
capable of building the intense spot in between the other 
two. However, in the coherent case £ = 1, shown in Fig. El 
(d)-(f), beams focus and overlap less, and the beams to 
the right (e) are expelled from the region between the 
other beams. Also, the time scale of build-up dynamics 
is shorter for coherent then for incoherent beams. 
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FIG. 5: Temporal evolution of the output intensity distribution of the two-lobe left-propagating beam at the left face of the 
crystal. Dashed line at t = 116to shows the slice corresponding to Fig. 01(c), where the modulational instability breaks the 
transverse symmetry. 



Of special interest are those self-trapped structures 
that dynamically do not converge to a steady state. Such 
structures represent novel time-dependent, as well as z- 
dependent waveguides that can not be described by the 
usual steady-state theory of spatial solitons. Whereas the 
^-dependence can be ascribed to the general definition of 
longitudinal waveguide modes, the time-dependence is a 
novel feature, caused by the slow response of PR crystals. 
An example is depicted in Fig. 0] where a collision of 
three against two power-matched coherent beams is pre- 
sented. The initial configuration is such that the three 
beams propagating to the right interfere constructively 
(b), to overlap with the two counterpropagating out-of- 
phase beams (c). These two beams, propagating to the 
left, repel and overlap with the two outside-lying oppo- 
site beams. During the time evolution of this dynami- 
cal state we have observed several alternations of trans- 
versely symmetrical structures, similar to the one shown 
in Fig. 01 and identified such behavior as a quasi-periodic 
self-oscillation Q, clearly seen in Fig. for t < 116to. 
At that point the development of transverse symmetry- 
breaking instability is observed, which results in irregular 
dynamics, shown in Fig. for t > 116to. 



In conclusion, we have developed a theory of self- 
trapped bidirectional waveguides. In counterpropagat- 
ing geometry, the inclusion of time dependent effects was 
found to be crucial for the formation of joint waveguiding 
structures. We demonstrated the generation of a counter- 
propagating (1+1)D vector soliton numerically and pro- 
posed a more general class of non-soliton steady-state 
solutions. The level of temporal coherence of interacting 
beams influences the mutual coupling due to the forma- 
tion of a refractive index grating. In addition to the gen- 
eration of steady-state induced waveguides, the dynamic 
alternation of states following by transverse modulational 
instability, as well as the onset of longitudinal modula- 
tional instability were observed. 
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